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ABSTRACT 



Context. The existence and origin of large spatial temperature fluctuations in H n regions and planetary nebulae are assumed to 
explain the differences between the heavy element abundances inferred from collisionally excited and recombination lines, although 
this interpretation remains significantly controversial. 

Aims. We investigate the spatial variation in electron temperature inside NGC 346, the brightest H n region in the Small Magellanic 
Cloud. 

Methods. Long slit spectrophotometric data of high signal-to-noise were employed to derive the electron temperature from measure- 
ments derived from localized observations of the [Om](/14959 + /15007)//14363 ratio in three directions across the nebula. 
Results. The electron temperature was estimated in 179 areas of 5"xl'.'5 of size distributed along three different declinations. A 
largely homogeneous temperature distribution was found with a mean temperature of 12 269 K and a dispersion of 6.1%. After cor- 
recting for pure measurements errors, a temperature fluctuation on the plane of the sky of f 2 = 0.0021 (corresponding to a dispersion 
of 4.5%) was obtained, which indicates a 3D temperature fluctuation parameter of t 1 » 0.008. A large scale gradient in temperature 
of the order of -5.7 ± 1.3 K arcsec -1 was found. 

Conclusions. The magnitude of the temperature fluctuations observed agrees with the large scale variations in temperature predicted 
by standard photoionization models, but is too small to explain the abundance discrepancy problem. However, the possible existence 
of small spatial scale temperature variations is not excluded. 

Key words. ISM: Hn regions - ISM: individual objects: NGC 346 



1. Introduction 

The most important unsolved issue in the study of photoionized 
nebulae is why the chemical abundances of heavy elements de- 
rived from recombination lines are systematically higher than 
those derived from collisionally excited lines, the so-called 
"abundance discrepancy" (AD) problem. For Hn regions, the 
ratio between the abundances e stimated from permitted and for - 
bidden lines is about ADF w 2 (Garcfa -Roias & Estebanll2007l) . 
Much larger abundance discrepa ncies (2 $ ADF < 2 0) have 
been found in planetary nebulae dLiu et al. L I2000U2001I) . Spatial 
temperature fluctuations have been proposed to explain the AD 
pr oblem. Th e se tem perature fluctuations were initially proposed 
bv lPeimber^ dl967l) to account for the considerable discrepancy 
found between the temperature estimates for Hn regions ob- 
tained by different methods. However, the high levels of temper- 
ature fluctuation required to solve the AD problem are not pre- 
dicted by standard photoionisation models ( Kingdon & Ferland, 
1 1995b . Chemical or density inhomogeneities of different charac- 
teristics and origins have been proposed to explain the high tem- 
perature fluctuations needed 



Kingdon & FerlandL 119981: 



2007: lTsamis et all 



d dGiammanco & Beckman. 2005; 
iLiuetall l2000t [ Stasinska et all 
94), 



2004; Viegas & Clegg, 1994), but the results 



of analysis based on observations are still inconclusive. 



Direct determinations of electron temperature fluctua- 
tions in the plane of the sky have been obtain e d for 
abo u t a do z en of planetary ne b ulae dKrabbe & CopettiL [2005; 
lEiul 119981; iRubin et all 120021; iTsamis et all l2008t) a nd two 
Hn regions, the Orion Nebula dMe sa-Delga do et all 120081 ; 
lO'Dell et all 120031; IRubin et all 120031) and the 30 Doradus 
Nebula dKrabbe & Copettil |2002|) . Localized measurements of 
electron temperature were obtained across the object surface 
area and small temperature fluctuations were inferred. 

In this paper, we describe a study of spatial variations in elec- 
tron temperature within NGC 346, the brightest H n region in the 
Small Magellanic Cloud (SMC). This object is also designated 
as LHA 115-N66 in the catalogue o f emission nebulae in the 
Magellanic Clouds by iHenizd (1 19561) and DEM S 103 in that of 
Davi es et all d!976l) . or more concisely as N66 and DEM 103, 
respectively. NGC 346 is a 14'xll' nebula with complex struc- 
ture, composed of arcs, filaments, and compact blobs of about 
30" or less (e.g. the component N 66A). The temperatures were 
derived from localized measurements of the [Om](/14959 + 
/15007)//14363 emission line ratio obtained from long-slit spec- 
trophotometric observations of high signal-to-noise ratio. 



2. Observations and data reduction 



Send offprint requests to: Vinicius de A. Oliveira, 
e-mail: vinic99@mail.ufsm.br 



The observations were performed on the nights of 8 and 9 
November 2002 at the Laboratorio Nacional de Astroffsica 
(LNA), Brazil, with the Cassegrain spectrograph attached to 
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Table 1. [O in] ratio and electron temperature statistics 
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Fig. 1. Spacial profile of the H/3 flux (in the units 
KT^ergscirrV 1 ), the [Om](^4659 + /15007)M4363 ra- 
tio (R), and the derived electron temperature along the A<5 = 
30"N. The position in horizontal axis is relative to the reference 
star. 



the 1.6 m telescope. The detector used was a SITe CCD of 
2048 x 2048 pixels. The spatial scale was 0756pixel _1 . Using 
a grid of 1200 grooves mrrr 1 , we obtained spectra covering 
the wavelength range of 4050 A to 5030 A with a dispersion 
of 0.5 A pixel 1 and resolution of 2.7 A. The slit used had 
an entrance of 360" x 1'.'5 on the plane of the sky and was 
aligned along the east-west direction. The slit was centered at 
3 different declinati ons relative to the reference star MPG 470 
dMassev et all fl989h (a = h 59 m 06 s and 6 = -72°10'34", 



J2000): A6 = 0", A6 = 30" N, and A6 = 30" S. At declination 
offset A6 = 30" N, the center of slit was moved almost 70" to the 
east of the reference star to cover a brighter section of NGC 346. 
To limit the effects of cosmic rays, three exposures of 20 minutes 
were completed at each slit position. At A5 = 30" S, two addi- 
tional exposures of 10 minutes were acquired. Several bias and 
dome fiat-fields exposures were taken at the beginning and the 
end of each night. For wavelength calibration, spectra of He-Ar 
lamp were taken before and after each object exposure. For flux 
calibration, the spectrophotometric standard stars HR9087 and 
HR 1544 were observed. 

Data reduction was performed using software, and 

the standard procedures for bias correction, flat-fielding, cosmic- 
ray event removal, and wavelength and flux calibration were 
followed. From the mean of the two-dimensional spectra ob- 
tained at each declination offset A6, we extracted a series of 
one-dimensional spectra from contiguous regions in length of 5" 
along the slit axis. To ensure that independent ID spectra were 
extracted from 2D spectra along the given north-south strips, we 
defined fiducial positions along the slit axis; this was achieved 
by identifying, inside the 2D spectra spatial profiles, the detected 
star position and comparing this with the coordinates of this star 
in direct images of the region. This procedure generated 179 one- 
dimensional spectra from regions of 5"xl'/5 of size. 

The line fluxes were measured by Gaussian fitting of the line 
profile over a linear local continuum. The measurements were 
completed with the splot routine of the IRAF package. The error 
associated with the line fluxes were estimated to be <x 2 = crL + 
cr^ ont , where <T\{ ne is the Poisson error of the emission line and 
cr cont is the error due to the continuum noise, calculated to be 
c C ont = yNAcTr.m .s ., where N is the number of pixels covered 
by the emission line, A is the dispersion of the spectrum (units 
of wavelength per pixel), and cr l m s is the root mean square of 
the continuum flux density (flux per unit wavelength). All line 
intensities were normalized to Hfi and corrected for the effect of 
interstellar extinction, by comparing the observed ratios Hy/H/3 
with the theoretical one, calculated by Storev & Hummer ( 1995) 
for an electron temperature of 10000 K and a de nsity of 100 
cm" 3 . The reddening law for the SMC obtained by iPrevot et alj 
d 19841) was used. 



3 IRAF is distributed by the National Optical Astronomy Observatory 
(NOAO), which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc., under cooperative agreement 
with the National Science Foundation. 
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Fig. 2. Same as Fig.Q] but for A6 = 0". 
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Fig. 3. Same as Fig. [I] but for AS = 30"S. 
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3. Determination of the electron temperature 

Electron temperatures were calculated from the [Om](/i4959 + 
/l5007)//i4363 line intensity ratio by solving numerically the 
equations of equilibrium for the n-level atom (« = 6) using the 
temden routine of the nebular package of STSDAS/IRAF (see 
IShaw&Dufourlll995h . The assumed values for energy levels, 
transition pro b abilities and collis ion s trengths were taken from 
lBoweddl96db . lWiese et al.l(ll996l) . and lLennon & Burkd(ll994l) . 
respectively. The dependence of the calculated electron temper- 
ature T e on the assumed electron density N e is very weak. In 
our estimates, we adopted a fixed electron density of 100 cm -3 
as a representative value. Electron de nsities in the range of 50 
to 500 cm' 3 were found for NGC 346 Jbufour & Hariowl 1 19771 
iPeimbert et all 120001: iTsamis et all 120031) . We calculated that, 
for electron densities in this range, the errors in the electron tem- 
perature estimates due to uncertainties in density would be below 
0.1%. 



4. Results and Discussions 

Figures QJ[3] show the spatial profiles along the slit of the mea- 
surement of H/3 flux, [Oin] ratio and electron temperature. 
Table Q]presents a summary statistics of the [Oin] ratio and the 
inferred electron temperature, including the number N of dif- 
ferent nebular sections (one-dimensional spectra), the minimum 
and maximum, the median, the first quartile Ql (upper limit of 
the 25% lowest values) and the third quartile Q3 (lower limit of 
the 25% highest values). We also indicate both the mean and 



the standard deviation cr weighted by the H/3 flux. The elec- 
tron temperature estimates are quite homogeneous. The scatter 
is small with 50% of the values less than 4.4% from the overall 
median temperature of 12 331 K. For a perfect Gaussian dis- 
tribution, the standard deviation is related to the quartiles by 
cr w 1.35(23 - Ql). Assuming that this relationship holds for 
the observed temperature distribution, we estimate a dispersion 
of 6.5% from the non-parametric statistics. Part of the scatter is 
real and part is due to measurement errors, as can be inferred by 
the increase in scatter at the fainter areas at the edges of the neb- 
ula. From the parametric statistics (mean and standard deviation 
weighted by the H/3 flux), which provided more weight to the 
more precise data, we estimated a weighted mean electron tem- 
perature of 12 269 K and a weighted standard deviation of 6. 1 %. 
The similarity between these two independent estimates of the 
dispersion in the data indicates that the measurement errors are 
not dominant sources of scatter in the temperature data. 

For comparison with our results, we recalculated the elec- 
tron temperature from line intensity ratios collected from the lit- 
erature. Table [2] presents the temperature diagnostics used, the 
line ratio values, the electron temperatures recalculated with up- 
date atomic parameters, the references for the observational data, 
and the authors's label for the extraction window for these data, 
which are shown in Fig. |4] The temperature estimates obtained 
from the data in the literature, although corresponding to differ- 
ent areas on the plane of the sky and different ionization zones, 
are very similar to both each other and our values. For example, 
the [O m] temperatures for data from the literature have a mean 
of 12685 K (3% higher than our mean value) and a dispersion 
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Table 2. Electron temperatures calculated using data from the 
literature 



Ratio 


r e (K) 


Ref./Region 


Ratio 


r e (K) 


Ref./Region 


[Ora] (,14959 + ,15007)/,14363 


99.8 


12 875 


[1] 


105.1 


12 625 


[2] 


107.4 


12519 


[3] 


110.9 


12 370 


[4] NW 


117.9 


12 091 


[4] SE 


96.6 


13 044 


[5] I 


96.7 


13 040 


[5] II 


90.1 


13 520 


[6] 1-lc 


92.6 


13 263 


[6] l-2b 


130.3 


11665 
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12 245 
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85.7 
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[7] 1 
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13 279 
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87.2 
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113.2 
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111.1 


12 363 


[7] 13 


108.0 


12492 


[7] 14 


105.5 


12 606 


[7] 15 


99.1 


12916 


[7] 16 


95.3 


13 113 


[7] 17 


89.1 


13 472 


[7] 18 


108.2 


12487 


[8] 








[Nn] (,16548 + /16583)//15755 


52.22 


13510 


[8] 








[On] (,13726 + ,13729)/(,17320 + ,17330) 


40.8 


11627 


[5] I 


35.5 


12 747 


[5] II 


43.0 


11707 


[7] 1 


45.6 


11305 


[7] 2 


35.4 


13 324 


[7] 3 


35.7 


13 262 


[7] 4 


45.8 


11263 


[7] 5 


51.8 


10 501 


[7] 11 


41.5 


11974 


[7] 12 


37.7 


12 768 


[7] 13 


34.0 


13731 


[7] 14 


24.9 


17 902 


[7] 15 


32.7 


14 125 


[7] 16 


42.6 


11773 


[7] 17 


33.7 


13813 


[7] 18 


36.7 


12459 


[8] 


[Sn] (,16716 + ,16731)/(,14068 + ,14076) 


15.2 


8613 


[7] 12 


15.5 


8 525 


[7] 13 


14.0 


9110 


[7] 17 


10.8 


10 792 


[8] 



HI lAller & Faulkner! dl962h: \2] 



Refere nces: 

dl964l): T31 iDufouri j!975l); [4] 
Peimbert & Torres-Peimbert] (|1976); 



Peimbert et alj d200d) : r81 lTsamis et ai1d26o3¥ 



Dufour & Harlow 
T61 |Pagel et alj 



(1977); 
G978I); 



[5] 
[7] 



of only 3.9%, while the mean [On] temperature is 12 153 K (1% 
lower than our mean value). A significant source of the discrep- 
ancies between different temperature estimates is error in the in- 
dividual measurements, which are at their lowest of the order of 
2-3%. 

All observations support the conclusion that NGC 346 has 
a homogeneous temperature structure with variations of about 
+6% or less. It is true that observations, which collect light 
that has been integrated along the line of sight, tend to smooth 
out small spatial-scale fluctuations of any line ratio, but local- 
ized measurements are able to resolve global internal gradi- 
ents. Any significant large-scale systematic variation in elec- 
tron temperature should then be revealed by high signal-to- 
no ise ratio observa tions, as in studies of the Orion Nebula 
alter et all. 1 19921) a nd th e planetary nebulae NGC 6720 
Garnett & DinersteinL l200ll) . NGC 4361 fl998l) . and 



NGC 2438, NGC 2440, NGC 3132, N GC 3242, NGC 6302, and 
NGC 7009 dKrabbe & CooettiL 120051) . In NGC 346, we found, 
along the declination offsets AS = 0" and A<5 = 30"N, a small 
but clear electron temperature gradients of -5.7 ± 1.3 K arcsec -1 , 
which is statistically sign ificant at the 95% co nfidence level. 
Using the data obtained by Peimb ert et alj (Hp00) for areas along 
a declination offset close to AS = 0" (their regions 11-18; see 
Fig.|H, we found a similar temperature gradient of -5.2 ± 0.8 K 




20 s 59 m Cf 
a (2000) 



Fig. 4. Apert ures of the obse rvations o f NGC 346 by: (OCR) thi s 
paper, ( D75) IDufouri d 19751). (DH77) iDufour & Harlowl d 19771) 
(PTP7 6) IPeimbert & Torres Peimberll dl976l) . (P7 8) iPaeel et at 
dl978|) . (POO) IPeimbert et alj d200oT) . and (T03) iTsamis et al 
320031) . 
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Fig. 5. Linear regression analysis for the data along declina- 
tion offset AS = 0"( squares, solid line) and for the data by 
IPeimbert et al] (120001) (triangles, dashed line). The curved lines 
define the 95% confidence bands for the regression lines. 



arcsec 1 . Figures [5] shows the result of a linear re gression analy- 
sis of our data and a comparison with the data of Peimbert et al. 
(|2000l) . Along the declination offset A5 = 30" S, no significant 
gradient was found. 



4.7. Magnitude of the electron temperature fluctuations 

The discrepancy between the temperature estimates obtained by 
means of forbidden line ratios and thos e derived from th e Balmer 
and radio continuum was explained bv lPeimbertl d 1 9671) as a con- 
sequence of spatial electron temperature fluctuations inside the 
nebula. To quantify the effects of these fluctuations on the esti- 
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mated temperature, Pei mbertl ( 1 1967b introduced the temperature 
fluctuation parameter f 2 defined to be 



f(T e -T ) 2 NiN e dV 

T 2 jNiN e dV 



with 



To 



J T e NiN e dV 
fNiN e dV 



(1) 



(2) 



where Ni is the density of the ion used to measure the temper- 
ature, integrations are calculated over the volume V of the neb- 
ula. Therefore, the parameters To, t 2 , and Vf 2 are the mean, the 
relative variance, and the relative standard deviation of the tem- 
perature distribution weighted by the square of the local density, 
respectively. The t 1 parameter has also been used to correct the 
chemical abundances obtained from collisionally excited lines 
for the effects of temperature fluctuations. Usually, ad hoc val- 
ues for f 2 are used (including t 2 = 0), or t 2 is estimated indirectly 
from the comparison between temperature values obtained from 
different indicators. 

A direct estimation of t 2 can be obtained from localized mea- 
surements of the electron temperature across the nebula. We may 
rewrite Eq.Q]and find a lower bound for t 2 as 



2 f <(r e - To) 2 } a EadQ. f «r e >n - T fE a dQ. 
t 2 = ^ > ^ ^ 4 r (3) 



T 2 J End®. 



Tq J EcidQ. 



where: (.)n represents the mean value weighted by the square 
density along the direction £2; E& = J NjN e dl is the emission 
measure; dQ. and dl are respectively the elements of solid angle 
and distance along the line of sight. Therefore, the f 2 ky parameter 
is the relative variance of the continuous distribution of (T e )a, 
the mean temperature along the line of sight at the direction Q 
weighted by the emission measure En- 



(1998), 



A discrete approximation for f 2 ky , first proposed by |Li 



r(obs) = z , 



(4) 



where T' e and Fj(HJ3) are the electron temperature and the H/3 
flux obtained for the aperture i, respectively. This quantity rep- 
resents the variance (relative to the square of the mean temper- 
ature) of the temperature values measured at different apertures, 
weighted by the H/? flux. Because part of this variance is due 
exclusively to errors in the measurements, the final estimation 
of t 2 should be corrected by t 2 = f 2 (obs) - t 2 noK , where f 2 Tors 
is the quadratic mean of the relative errors of temperature, also 
weighted by the Hf3 flux. 

Since the temperature measured for any aperture is the 
mean value inside a volume crossing the nebula, any temper- 
ature variations within distance scales smaller than the aper- 
ture size should be smoothed out. Hence, t\ is a biased esti- 
mator of f 2 , such that t\ < f 2 ky < t 2 . However, localized fluc- 
tuations in temperatures are unexpected (see Ferlandl 12001 1) . 
Photoionization models of chemically homogeneous nebulae in- 
ste ad predic t smooth temperature gradients. 

ICopettil d2006l) demonstrated with simulation that localized 
temperature estimates obtained for ions occurring in significant 
parts of the nebula, such as the [O in] and the Balmer jump tem- 
peratures, are good tracers of the internal gradient of tempera- 
ture, and that the values of t\ calculated from these temperatures 



contribute significantly to the total variance, with f 2 (Om)/f 2 ~ 
25% for a typical H n region. 

For NGC 346, we obtained an observed temperature vari- 
ance of ? 2 (obs) = 0.00375 (corresponding to a dispersion of 
6.1%), a mean quadratic error of 4.1%, and a corrected vari- 
ance in the temperature distribution on the plane of the sky of 
f 2 = 0.0021 (or equivalently, a dispersion of 4.5%). This value 
is similar to t 2 - 0.0025 obtained for the 30 Doradus Nebula 
dKrabbe & Copettil |2002|) and is i nside the range of 0.001 1 to 
0.0050 found in planetary nebulae dKrabbe & Copettil 120051) . 

Using the [O ml U 4959 + A5001) I AA2>62> ratios measured by 
IPeimbert et ail d2000l) for 13 areas of NGC 346 and recalculat- 
ing the electron temperature with the same atomic data used 
in the present paper (see Table [2J, we calculated a variance 
in the temperature distribution on the plane of the sky of only 
f 2 = 0.00062, which is even lower than the value that we esti- 
mated from our data, possibly because their data extraction win- 
dows were up to a factor of 7 larger than ours and were located 
within a smaller area of NGC 346. 

Us ing the approximation f 2 (Oin)/f 2 ~ 25% from ICopettil 
d2006l) . we estimate a global temperature fluctuation parameter 
of t 2 w 0.008 for NGC 346, which is entirely compatible with 
the values predicted by photoionization models for an object 
with the characteristics of NGC 346, b ut is too small to e xplain 
the abundance discrepancy problem. Tsam iset al.l d2003l) deter- 
mined the ++ /H + abundance ratio for NGC 346 from collision- 
ally excited and recombination lines and found an abundance 
discrepancy factor of ADF = 2.3; this would require signifi- 
cant temperature fluctuations, corresponding to f 2 = 0.09, to be 
described by the temperature fluctuation scenario. Although we 
cannot exclude the possibility of small spatial-scale temperature 
fluctuations, since they may not be detected by our localized ob- 
servations, it is difficult to accept that variations in temperature 
of dispersion around the mean of the order of 30% would remain 
undetected. 

We emphasize that although the parameter f 2 (0 in) probes 
the temperature variation in the ++ zone, the correction applied 
to calculate the 3D temperature fluctuation parameter t 2 takes 
into account both the projection effect of the temperature mea- 
surements on the plane of the sky and the variation in temper- 
ature outside the ++ zone. This correction is based on model 
simulations, and is therefore model dependent. In particular, for 
a s pherical symme tric and constant density nebula, the models 
by ICopettil (120061) predict values of f 2 (Om)/f 2 of between 8 
and 40% for a wide range of the input parameters (hydrogen 
density A^h, effective temperature of the ionizing star T e s, ion- 
ization parameter U). For a more restricted range of parameters 
(40000$ 7/ eff £50 000 K,N H ~ 100- 1000 cm" 3 , log U * -3), 
corresponding to typical H n regions like NGC 346, a more lim- 
ited range of ratios of 15% < f 2 (Om)/f 2 < 35% is found. We 
therefore expect an error of about 40% for the estimate of t 2 
obtained from f 2 (0 in), due solely to the model dependency on 
the input parameters. Therefore, this estimation of the temper- 
ature fluctuation parameter f 2 , and any other obtained by any 
method, should always be considered as approximate estimation 
because of the large associated errors. Even though, it is clear 
that the temperature fluctuation parameter of f 2 = 0.09 obtained 
from the comparison between abundances derived from recom- 
bination and collisionally excited emission-lines is incompatible 
with the order of magnitude smaller value of t 2 = 0.008 obtained 
from localized determinations of electron temperature across the 
ne bula. 

IPeimbert et all (120001) inferred a value of f 2 = 0.022 ±0.012 
for NGC 346 (and a minimum of t 2 = 0.0013) by comparing the 
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mean temperatures measured using different indicators, which 
estimates the temperature fluctuation that accounts directly for 
the temperatures differences between different ionization zones. 
This value is larger but compatible within the error bars with 
our own estimation. It is still too small however to explain the 
abundance discrepancy problem. 

5. Conclusions 

We have used long-slit spectrophotometric data of high signal- 
to-noise ratio in the range of 4050 - 5030 A to study the spa- 
tial variation in electron temperature in the H n region NGC 346. 
Electron temperatures were derived from the [Om](i4959 + 
^5007)/i4363 line ratio for 179 areas of 5" x 175 distributed in 
three different declination. Our main results are the following: 

1. A largely homogeneous electron temperature distribution 
was found with a mean value of 12 269 K and a standard 
deviation of 6. 1 %, both weighted by the H/? flux. 

2. A small, but statistically significant, large scale gradient in 
temperature of the order of -5.7 ± 1.3 K arcsec -1 was found 
along the east-west direction at two different offset positions. 

3. About half of the variance of the temperature was attributed 
to measurement errors. The remainder corresponded to a 
temperature fluctuation on the plane of the sky of t\ = 0.0021 
(equivalent to a dispersion of 4.5%), which implied a 3D 
temperature fluctuation parameter of f 2 « 0.008. 

4. The magnitude of the temperature fluctuations observed is 
in agreement with the large scale variations in temperature 
predicted by standard photoionization models, but is too low 
to explain the abundance discrepancy problem. However, the 
existence of small spatial scale temperature variations is not 
excluded by the present observations. 
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